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ABSTRACT 

At redshift z > 5 Type II supernovae (SNII) are the only known dust sources with 
evolutionary timescales shorter than the Hubble time. We extend the model of dust 
formation in the ejecta of SNII by Todini & Ferrara (2001) to investigate the same 
process in pair-instability supernovae (PISN), which are though to arise from the 
explosion of the first, metal free, very massive (140-260 Mq) cosmic stars. 

We find that 15%-30% of the PISN progenitor mass is converted into dust, a value 
> 10 times higher than for SNII; PISN dust depletion factors (fraction of produced 
metals locked into dust grains) range between 0.3 and 0.7. These conclusions depend 
very weakly on the mass of the PISN stellar progenitor, which instead affects consid- 
erably the composition and size distribution. For the assumed temperature evolution, 
grain condensation starts 150-200 days after the explosion; the dominant compounds 
for all progenitor masses are Si02 and Mg2Si04 while the contribution of amorphous 
carbon and magnetite grains grows with progenitor mass; typical grain sizes range 
between 10"'^ and a few O.l/im and are always smaller than l/xm. We give a brief dis- 
cussion of the implications of dust formation for the IMF evolution of the first stars, 
cosmic reionization and the intergalactic medium. 
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In the recent years, dust has been recognized to have an in- 
creasingly important role in our understanding of the near 
and distant Universe. The dramatic effect of dust at low and 
moderate redshifts has been immediately recognized when a 
recostruction of the cosmic star formation history from rest- 
frame UV/ visible emission was first attempted: dust grains 
absorb stellar light and re-emit it in the FIR. Thus, even 
a tiny amount of dust extinction can lead to a severe un- 
derestimate of the actual star formation rate (Pettini et 
al. 1998; Steidel et al. 1999). New IR, FIR and submm fa- 
cilities have revealed the existence of populations of sources, 
such as SCUBA « > 1 sources, that are thought to be 
dust-enshrouded star forming galaxies or AGNs (Small et 
al. 1997; Hughes et al. 1998) or the Extremely Red Objects, 
which are at least partly populated by dusty star-forming 
systems at 2: ~ 1 (Cimatti et al. 2002). Finally, dust plays 
a critical role in galaxy evolution, accelerating the forma- 
tion of molecular hydrogen (H2), dominating the heating of 
gas through emission of photoelectrons in regions where UV 
fields are present and contributing to gas cooling through 
IR emission (see Draine 2003 for a recent thorough review) . 

These observational evidences have motivated a series 



of studies aimed at including a treatment of dust forma- 
tion within galaxy evolution models (Granato et al. 2000; 
Hirashita & Ferrara 2002; Morgan & Edmunds 2003). Since 
dust formation in the interstellar medium (ISM) is extremely 
inefficient (Tielens 1998), the preferential sites of formation 
are considered to be the atmospheres of evolved low-mass 
(M < 8Mq) stars, from where it is transported into the ISM 
through stellar winds (Whittet 1992). However, this mech- 
anism can not explain the presence of dust at redshifts > 5 
because at these high redshifts the evolutionary timescales 
of low-mass stars (0.1-1 Gyr) start to be comparable with 
the age of the Universe. 

Evidences for the presence of dust at high redshifts 
come from observations of damped Lya systems (Pettini 
et al. 1994; Prochaska & Wolfe 2002; Ledoux, Bergeron & 
Petitjean 2002) and from the detection of dust thermal emis- 
sion from high redshift QSOs selected from the SDSS sur- 
vey out to redshift 5.5 and re-observed at mm wavelengths 
(Omont et al. 2001; CariUi et al. 2001; Bertoldi & Cox 
2002). Very recently, Bertoldi et al. (2003) have reported 
the observations of three z > 6 SDSS QSOs at 1.2 mm, de- 
tecting thermal dust radiation. From the IR luminosities, the 
estimated dust masses are huge (> 10* M©) implying a high 
abundance of heavy elements and dust at redshifts as high 
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as 6.4 that can not be accounted by low-mass stars. Thus, 
dust enrichment must have occurred primarily on consider- 
ably shorter timescales in the ejecta of supernova explosions 
(Dwek & Scalo 1980; Kozasa, Hasegawa & Nomoto 1989; 
Todini & Ferrara 2001; Nozawa ct al. 2003). 

The strongest evidence for dust formation in supernova 
explosions was seen in SN 1987A (Moseley et al. 1989; 
Wooden 1997). In this event, the increased IR emission was 
accompanied by a a corresponding decrease in the optical 
emission and the emission-line profiles were observed to shift 
toward the blue (McCray 1993 and references therein). In 
another supernova explosion, SN 1998S, the observed evolu- 
tion of the hydrogen and helium line profiles argues in favor 
of dust formation within the cjocta as the redshifted side 
of the profile steadily faded while the blueshifted side re- 
mained constant (Gerardy et al. 2000). Dust emission has 
been seen in the supernova remnant Cassiopea A: similarly 
to SN 1987A, the total dust mass derived from IR lumi- 
nosities is less than expected from theory, suggesting that a 
colder population of dust grains may be present that emit 
at longer wavelengths and it is not detectable in the IR 
(Morgan & Edmunds 2003; Dunne et al. 2003). Finally, in 
a considerable rmmbcr of cases, supernovae have shown IR 
emission that was stronger toward longer wavelengths (see 
Gerardy et al. 2002 and references therein) . This IR excess 
has been generally interpreted as due to thermal emission 
from dust forming in the ejecta but alternative explanations 
exist; in particular, new IR observations of five Type II SNc 
have shown late-time emission that remains bright many 
years after the maximum and that it is hard to reconcile 
with emission from newly formed dust; IR echos from pre- 
existing dust in the circumstellar medium heated by the su- 
pernova flash might represent an alternative interpretation 
(Gerardy et al. 2002). 

Theoretical studies have started to investigate the pro- 
cess of dust formation in expanding SN ejecta. Most of the 
available models are based on classical nucleation theory and 
grain growth (Kozasa & Hasegawa 1987; Todini & Ferrara 
2001) with the exception of the computations recently per- 
formed by Nozawa et al. 2003, who have made an extensive 
analysis of dust grain properties forming in zcro-metallicity 
SNe of various masses. The model developed by Todini & 
Ferrara (2001) is able to predict the dust mass and proper- 
ties as a function of the initial stellar progenitor mass and 
metallicity. In spite of the many uncertainties and approx- 
imations, this model has been shown to satisfactorily re- 
produce the observed properties of SN 1987A and of the 
young dwarf galaxy SBS 0335-052 (Hirashita, Hunt & Fer- 
rara 2002). 

In this paper, we investigate the formation of dust 
grains in the ejecta of very massive primordial supernovae, 
which are commonly known as pair-creation or pair insta- 
bility SNe (PISN). These violent explosions are thought to 
terminate the life of metal-free stars with initial mass in the 
range UQMq < M < 260Mq (Heger & Woosley 2002). In- 
deed, detailed theoretical modelling of the nucleosynthesis 
and internal structure of these very massive stars has shown 
that after central helium burning electron/positron pairs are 
created, converting a large fraction of internal energy into 
rest mass of the pairs; as a consequence, the stars rapidly 
contract until explosive oxygen and silicon burning is able 
to revert the collapse and the stars are completely disrupted 



in giant explosions (Fryer, Woosley & Heger 2001). Below 
and above the mass range of PISN, the most likely outcome 
of the evolution of metal-free stars appears to be black hole 
formation, cleanly separating the contribution to metal and 
eventually dust production of such very massive stars from 
that of other mass ranges. 

We base our analysis on the model developed by To- 
dini & Ferrara (2001) and apply it to PISN. We therefore 
change the initial chemical composition of the ejecta and 
also the tliermodynamical/dynamical properties which de- 
termine the ejecta evolution. A larger amount of dust is 
expected to form in the ejecta of PISN with respect to Type 
II SNc, because of the larger amount of metals released and 
the absence of fallback of material onto the compact rem- 
nant (no remnant is expected to survive PISN explosions). 

The motivation for the present study comes from the 
increasing number of evidences which seem to indicate that 
the first stars that were able to form in the early Uni- 
verse from the collapse of metal-free gas clouds were indeed 
very massive, with characteristic masses of a few 100 Mq. 
These include detailed numerical simulations (Abel, Bryan 
& Norman 2000, 2002; Bromm, Coppi & Larson 1999, 2002; 
Brormn ct al. 2001: Ripamonti et al. 2002), semi-analytic 
models (Omukai & Nislii 1998; Nakamura & Umemura 2001, 
2002; Schneider et al. 2002) and several pieces of obser- 
vations (Hernandez & Ferrara 2001; Oh et al. 2001; Sal- 
vaterra & Ferrara 2002; Schneider et al. 2003). In particu- 
lar, some of these studies have pointed out that because of 
the reduced gas coohng efHciency, low-mass star formation 
is strongly inhibited before a minimum level of metal enrich- 
ment of the collapsing gas cloud has been reached (Bromm 
et al. 2001, Schneider et al. 2002). The value of this min- 
imum level, Zci is very uncertain, but likely to be between 
lO"*^ and IQ-^Zq (Schneider et al. 2002; 2004). Within 
this critical range of metallicities, the presence of dust ap- 
pears to have a major role, providing an additional pathway 
for cooling the gas, that fragments into lower mass clumps 
enabling the formation of second-generation low-mass stars 
(Schneider et al. 2003). Thus, to estimate the efficiency of 
this dust-regulated low-mass star formation channel in very 
metal-poor gas clouds, it is crucial to develop a model which 
is able to predict the amount of dust formed in the ejecta of 
PISN. 



2 DUST FORMATION MODEL 

The process of dust formation and the resulting dust grain 

properties depend on the physical conditions at the site 
of formation. In particular, several studies (see Kozasa, 
Hasegawa & Nomoto 1989 and references therein) have 
shown that the time of onset of grain formation depends on 
the temperature structure in the supernova ejecta whereas 
the grain composition mainly reflects its chemical composi- 
tion, which depends on the nucleosynthesis occurring during 
stellar lifetime {i.e. on the progenitor mass) and explosion. 

Thus, models of dust formation in supernova ejecta are 
based on specific prescriptions for the chemical composition 
and thermodynamics of the expanding gas. At the onset of 
shock generation, which occurs at the boundary of the in- 
nermost Fe-Ni core, the progenitor star is characterized by 
the standard stratified (onion-skin) structure. During shock 
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propagation through the star, the gas undergoes a new phase 
of (explosive) nucleosynthesis and mixing of the internal lay- 
ers is thought to occur at least up to the outer edge of the he- 
lium layer, as suggested by observations of early emergence 
of 7 and X-rays in SN 1987A (Kumagai et al. 1988 and 
references therein). These observations confirm that macro- 
scopic mixing of the ejecta occurs through Rayleigh- Taylor 
instabilities but it is still very debated whether mixing is 
extended at molecular level (Clayton et al. 1999; 2001). 

The model of Todini & Ferrara (2001) (but see also 
Kozasa et al. 1989) is based on the assumption that ma- 
terials are uniformly mixed from the center to the helium 
outer edge. As a consequence, the temperature and density 
within this metal-rich volume are assumed to be constant. 
When the shock reaches the surface of the progenitor, the 
star starts to expand and the expansion becomes homolo- 
gous. Thus, the velocity of gas at a given layer is constant 
in time and proportional to the radius from the center, i.e. 
V — R/t, where t is the time from the explosion and the 
expansion velocity is assumed to be 

where iJkin and Mtot are the kinetic energy and total mass 
ejected by the supernova. The temperature of the expanding 
ejecta is determined by various heating and cooling mech- 
anisms. Todini & Ferrara (2001) following Kozasa et al. 
(1989) assumed that when the gas reaches the photosphere, 
the gas temperature is equal to the photospheric tempera- 
ture and thereafter the temperature evolution follows from 
the assumptions of adiabatic expansion for a perfect-gas, so 
that 

/ y \ 3(1-7) 

r = T.(i + -t) . (2) 

where 7 is the adiabatic index and the quantities Ti and 
Ri are the photospheric temperature and radius obtained 
from the observational results of SN 1987A (Catchpole et 
al. 1987). The lack of observational constraints for PISN 
forces us to model the thermal evolution of the expand- 
ing ejecta on the basis of theoretical considerations. In par- 
ticular, the evolution and internal structure of PISN have 
been studied in great detail through numerical simulations 
(Fryer, Woosley & Heger 2001; Heger & Woosley 2002). 
Zero-metallicity progenitor stars with masses between 140 
Mq and 26OM0 evolve without significant mass loss until 
central helium burning. These stars, after central helium de- 
pletion, have enough central entropy that they enter a tem- 
perature and density regime in which electron/positron pairs 
are created in abundance, converting internal gas energy 
into rest mass of the pairs, without contributing much to 
the pressure. When this instability is encountered, the stars 
contract rapidly until explosive oxygen and silicon burning 
produce enough energy to revert the collapse and the stars 
are completely disrupted in a giant explosion. The maxi- 
mum central temperature and density during the bounce as 
a function of the PISN progenitor mass together with the 
elemental composition of the emerging ejecta are shown in 
Fig. 1. The dominant metal yields show that nucleosynthe- 
sis in pair-creation supernovae produce a total mass of O, 
Si, Mg and Al which is roughly independent of the initial 
progenitor mass but an Fe mass which varies greatly with 



m 
CO 

S 1 



: 1 ' ' 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




r .f^ 






^ — ' — * ^^^^ 


' • - 


- 




- 




/ Al 10= 




; / 


Fe / 


: 










1 1 1 1 1 1 1 1 1 1 1 1 1 1 






















" ' u 


V 




n r / 1 


. 1 . . . 1 . . . 1 . . . 1 . < < 1 . 


, , 1 T 


140 


160 180 200 220 240 
Progenitor Mass [Mg] 


260 



Figure 1. Top panel: dominant metal yields [Mq] of PISN 
(Heger & Woosley 2002) as a function of the progenitor stellar 
mass, after all unstable isotopes have decayed. Bottom panel: ki- 
netic energy (in units 10^^ erg), velocity of the ejecta (in units 
10'^ km/s), central temperature (in units 10® K) and density (in 
units 10^ gr/cm^) as a function of the progenitor stellar mass (see 
text) . 



the mass of the progenitor, being almost negligible for initial 
stellar masses < 2OOM0. It is important to stress that this 
should be interpreted as the iron mass only after all unstable 
isotopes have decayed. Indeed, unlike Type II SN progeni- 
tors, PISN progenitors do not build up an iron core before 
the explosion and the final Fe mass is generated through 
the decay of ^^Co. We will consider this process in detail as 
this decay provides the relevant destruction process of CO 
molecules. 

The internal structure of the star at this maximum cen- 
tral density appears to have an approximately constant den- 
sity, pc and temperature, Tc, up to the outer edge of the 
helium core (Fryer et al. 2001). We can adopt these values 
as the initial density and temperature of the volume where 
the metals are uniformly mixed at shock generation. Indeed, 
the time it takes the shocks to propagate through the low- 
density hydrogen envelope and reach the surface of the pro- 
genitor, starting the stellar expansion, is negligible. Adopt- 
ing this normalization in Eq. 2, we can follow the evolution 
of the supernova ejecta. The results are shown in Fig. 2. In 
the top and medium panels we show the evolution of radius 
and number density of the ejecta of a 195 Mq pair-creation 
supernova as a function of the time elapsed from the ex- 
plosion. In the bottom panel, we show the corresponding 
evolution of the temperature. Two different regimes can be 
identified: in the initial phase of the evolution, the ejecta 
is radiation-dominated and we assume an adiabatic index 
7 = 4/3. This appears to be consistent with the internal 
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Figure 2. Dynamics and thermodynamics of the ejecta for a 
195 Mq PISN. The top, medium and bottom panels show the 
time evolution of the radius, Log [R/km], number density, Log 
[n/cm^], and temperature, Log [T/K], of the expanding gas. The 
temperatures corresponding to the onset of grain condensation 
(T < 2000 K) and to the final state followed by the model (T = 
500 K) are indicated with two horizontal lines. 



Figure 3. Time evolution of the dust mass in different solid com- 
pounds synthctized by a 149Mq {top panel), 195M0 (medium 
panel) and 250Mq {bottom panel) PISN. The different lines cor- 
respond to AC (dot-dashed), AI2O3 (long-dashed), Mg2Si04 
(dashed), Fe304 (solid) and Si02 (dot-long dashed). The dot- 
ted lines indicate the corresponding values of the temperature of 
the ejecta. 



structure of the 250 Mq star at maximum central density in 
the simulation of Fryer et al. (2001), which shows a density 
and temperature profile in the outer H-envelope compati- 
ble with 7 ~ 4/3. However, for t > 111.26 days (where 
111.26 days is ^^Co e-folding time) energy deposition from 
decaying radioactive elements is negligible and the gas cools 
adiabatically with 7 — 5/3. A comparable thermal history is 
found by Nozawa et al. (2003) solving the radiative transfer 
equation and taking into account the energy deposition by 
radioactive elements. For the same PISN progenitor mass, 
however, their normalization is slightly higher (by a factor 
~ 3) because of the different input stellar models (Umeda 
& Nomoto 2002). 



3 MODEL RESULTS 

In the model of Todini & Ferrara (2001), the formation of 
dust grains from a gaseous metal-rich medium is described 
as a two step process: (i) the formation of "critical clus- 
ters" at the corresponding condensation barrier and (ii)the 
subsequent growth of these clusters into macroscopic dust 
grains through mass accretion. A thorough description of 
the model, which follows the classic theory of nucleation, 
can be found in the paper of Todini & Ferrara (2001) to 
which we refer interested readers. 

In the present study, we follow the formation of seven 
solid compounds, namely, AI2O3 (corundum), iron, Fe304 



(magnetite), MgSiOs (enstatite), Mg2Si04 (forsterite) , Si02 
and amorphous carbon (AC) grains. The corresponding 
chemical reactions are listed in Table 1. Numerical constants 
can be found in Table 1 of Todini & Ferrara (2001) except for 
Si02, for whom we have adopted the values used by Nozawa 
et al. (2003). 

It is important to note that, as we will see, AC grains 
can form also if the ejecta composition is richer in oxygen 
than in carbon (O > C) as Clayton, Liu & Dalgarno (1999) 
have shown. The formation of CO and SiO molecules in SN 
ejecta can be very important for dust formation, because 
carbon atoms bound in CO molecules are not available to 
form AC grains and SiO molecules take part in the reactions 
which lead to the formation of MgSiOs, Mg2Si04 and Si02 
(see Table 1). Thus, the process of molecule formation in the 
expanding ejecta is followed at temperatures T < 2 x 10* K 
together with the ^''Co ^^'Fe radioactive decay as the 
impact with energetic electrons produced during this decay 
represents the main destruction process of CO molecules. 

The formation of dust grains is followed until the tem- 
perature of the ejecta has decreased to about 500 K, because 
at this stage condensation processes are terminated for all 
solid compounds. 

The time evolution of the dust mass synthetized in dif- 
ferent compounds is shown in Fig. 3 for three representa- 
tive PISN progenitor masses, 149M0, 195Af0 and 250Mq. 
Dust grains start to condensate when the temperature of 
the expanding gas has dropped below 2000 K, about 150- 
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Table 1. Chemical reactions included in dust formation calculations. 



Solid compound 


Chemical reaction 


AI2O3 


2A1+30 AI2O3 


Fe 


Fe(g) ^ Fe(s) 


Fe304 


3Fe+40 Fes04 


MgSiOs 


Mg+SiO+20 ^ MgSiOa 


Mg2Si04 


2Mg+SiO+30 Mg2Si04 


Si02 


SiO+O ^ Si02 


AC 


C(g) -> C(s) 



200 days after the explosion (see dotted lines). During this 
time interval, a fraction of the ^®Co initially present in the 
ejecta has decayed to Fe. Thus, at the onset of grain con- 
densation, some C atoms are bound in CO molecules and 
others are available to form amorphous carbon grains. The 
resulting dust mass is mainly composed by silicates and mag- 
netite grains. Amorphous carbon grains are the first to form, 
followed by corundum, forsterite, magnetite and ultimately 
Si02. This evolutionary sequence simply reflects the differ- 
ent condensation temperatures. The relative abundances of 
the different grain species depend on the progenitor mass: 
the I49M0 star ejects less iron (®^Co) than higher mass 
PISN and thus its final dust composition is characterized by 
a lower concentration of magnetite and carbon grains (the 
latter because most of the initial carbon grains are locked in 
CO molecules). 

3.1 Final dust masses 

The total dust mass formed in different solid compounds and 
the total mass of CO molecules synthetized in the ejecta as 
a function of the progenitor stellar mass are shown in Fig. 4. 
For all but the smallest progenitor star, silicate grains ap- 
pear to be the dominant dust compounds. Indeed, all SiO 
molecules initially present in the ejecta take part in the reac- 
tions leading to the formation of forsterite and Si02 grains 
and only CO molecules are left at the end of dust condensa- 
tion. At the lowest mass bin, I4OM0, the initial ''^Co mass 
is not large enough to favour the formation of magnetite 
and amorphous carbon grains and only AI2O3 and Mg2Si04 
grains are formed. At higher initial stellar masses, the larger 
^^Co mass ejected by the PISN favours the formation of car- 
bon and magnetite grains and the final mass of CO molecules 
is correspondingly reduced. Note that MgSi03 grains are 
never formed because of the highest condensation tempera- 
ture and largest nucleation current of Mg2Si04 grains, which 
lock all the Mg initially present in the ejecta. Similarly, Fe(s) 
are never formed because all the Fe which is produced by 
radioactive decays is locked into magnetite grains. Note that 
we do not consider the formation of MgO grains which, how- 
ever, can occur only in the ejecta of the 14OM0 PISN where 
Si is less abundant than Mg, and only 40% of the initial Mg 
mass is depleted onto Mg2Si04 grains. 

3.2 Grain size distribution 

The final grain size distributions are shown in Fig. 5 for 149 
Mq, 195 Mq and 250 M© PISN. The characteristic grain 
sizes range between 10~^ and a few 0.1 /im, depending on 
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Figure 4. Top panel: final dust mass of different solid compounds 
formed in PISN ejecta as a function of the mass of the stellar pro- 
genitors. Bottom panel: total mass of CO molecules synthetized 
in the ejecta as a function of the initial mass of the stellar pro- 
genitors. 



the grain species and on the mass of the progenitor star. 
Each dust grain grows by accretion to a final size which 
depends on the density regime when condensation occurs 
(higher densities favoring higher accretion rates) and on the 
abundance of the key species. Indeed, being the first to con- 
dense, amorphous carbon grains tend to have characteristic 
sizes larger than silicates and magnetite grains. AI2O3 grains 
are typically small because of the rather low-abundance of 
Al in the ejecta. Larger mass progenitor systems tend to 
have larger magnetite grains and smaller amorphous carbon 
grains. However, the properties of grain size distributions 
for different progenitor masses depend on the different ther- 
modynamics of the expanding ejecta which affects the rel- 
ative amplitude of cooling and nucleation rates. The model 
predicts that, for all progenitor masses, PISN lead to the 
formation of grains with radii smaller than 1 fim but sys- 
tematically larger than the typical sizes of grains formed 
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Figure 5. Final grain size distributions for different solid com- 
pounds synthetized by a 1A9Mq {top panel), 195M0 {medium 
panel) and 25OM0 {bottom panel) PISN. The distribution is com- 
puted when the temperature of the ejecta has decreased to 500 K 
and formation and accretion are terminated for all compounds. 



Figure 6. Top panel: total dust depletion factor defined as the 
dust-to-metal mass ratio as a function of the initial progenitor 
mass. Bottom panel: corresponding values of the ratio between 
the total dust mass synthetized and the initial progenitor mass. 



in SNII explosions (Todini & Ferrara 2001) with a reduced 
scatter among different grain species. 



3.3 Depletion factors 

The cosmological relevance of dust synthetized in PISN ex- 
plosions will depend mainly on the global properties of dust 
rather than on the nature and size of different compounds. 
In Fig. 6 we show the total dust depletion factor, defined 
as the dust-to-metal mass ratio released in the explosions, 
ct as a function of PISN progenitor mass. 
With the exception of the smallest mass bin, /dop ranges 
between 0.3 and 0.7 and tends to increase with increasing 
progenitor mass. In the bottom panel of the same figure, we 
show the ratio of the total dust mass and the initial pro- 
genitor stellar mass. For all but the smallest mass bin, the 
amount of dust formed is between 15 % and 30 % of the 
initial stellar mass. The peculiar behaviour of the 140 Mq 
progenitor is due to the original elemental composition of 
its ejecta, which is highly underabundant in Fe and Si and 
therefore only corundum is synthetized in the final ejecta 
(see also Fig. 4). 

For comparison, in Fig. 7 we show the same quantities 
but for SNII. Following Todini & Ferrara (2001), we show 
two cases corresponding to two different values of the total 
kinetic energy released in the explosion: case A corresponds 
to the low kinetic energy model (1.2 x lO^'^erg) and case B 
refers to the high kinetic energy model (1.9 x lO^^erg) (see 
also the original paper of Woosley & Weaver 1995) . We also 
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Figure 7. Same as in Fig. 7 but for Type II SNe. Case A and 
B refer to low (1.2 X lO^^erg) and high (1.9 X lO^^erg) kinetic 
energy scenarios and the dots (squares) represent the dust mass 
formed from Type II SNe with initial zero (solar) metallicity. 
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assume two different values for the initial metallicity of the 
progenitor stars, Z = (dots) and Z = Zq (squares). 

The moderate metal yields in SNII and the effect of 
fallback of material after the explosion onto the compact 
remnant lead to depletion factors which arc comparable to 
PISN, with values ranging from 20 % up to 70 % for Z = 
progenitors with masses < 22Mq. For larger masses, the 
depletion factor decreases in case B scenarios (higher kinetic 
energy) because of the larger amount of metals released. 
For case A scenarios (lower kinetic energy), instead, 25 Mq 
and SOM© SNII are predicted to have an /dep = 1- This is 
due to the fact that because of fallback, these stars eject 
only a rather small amount of metals, mostly in the form 
of carbon, which is completely depleted into AC grains. If 
the stellar progenitors have solar metallicities, the depletion 
factor ranges between 10% and 50% with a much reduced 
scatter between case A and B. 

In spite of these large depletion factors, the total mass 
of dust synthctizcd by PISN is significantly higher than that 
produced by SNII. Indeed, as already discussed in Todini & 
Ferrara (2001), in case A scenarios, Z = SNII synthetize 
a total dust mass which corresponds to a fraction between 
0.3 and 3% of the original stellar progenitor mass (thus, ^ 
0.1 — O.6M0 of dust per SN). These values are slightly larger 
if case B models are considered. When the initial stellar 
progenitors have solar metallicity, the resulting dust mass is 
typically a factor ~ 3 larger than for the metal-free case but 
it is always less than 6 % of the initial stellar mass. 



4 SUMMARY AND DISCUSSION 

We have investigated the process of dust formation in the 
ejecta of first stellar explosions, assuming that the first stars 
form with characteristic masses of a few 100 Mq and explode 
as pair-creation SNe. The study is based on an extension of 
the model developed by Todini & Ferrara (2001) for Type II 
SNe, which accounts for the different initial chemical compo- 
sitions of the ejecta and dynamical/thermodynamical prop- 
erties of the ejecta evolution. 

The main results of our analysis can be summarized as 
follows: 

(i) During pair-creation SN explosions, a significant 
amount of dust is synthetized out of the heavy elements 
produced in the progenitor stellar interiors during the main 
sequence lifetime and at the onset of the explosion due to 
explosive nucleosynthetic processes. The fraction of the orig- 
inal stellar mass which is converted into dust depends on the 
progenitor mass, with typical values ranging between 15% 
and 30%, resulting in SOA^q - GOA^q of dust produced per 
SN. These values are much larger than those found for Z = 
Type II SNe resulting from stars with masses between 12 
Mq and 40 Mq (O.IM© - O.6M0), even if these stars are 
assumed to be of solar metallicity (O.IMq - I.8M0). 

(ii) Because of the large amount of metals released in 
PISN explosions, dust depletion factors, defined as the dust- 
to-metal mass ratios, range between 30% and 70% depend- 
ing on the progenitor mass. 

(iii) The composition of the dust compounds depends 
critically on the thermodynamics and initial composition of 
the ejecta, i.e. on the progenitor mass. For the assumed 
temperature evolution, grain condensation starts 150-200 



days after the explosion and silicates, magnetite and AC 
grains are formed. The dominant compounds for all progen- 
itor masses are Si02 and Mg2Si04 while the contribution of 
amorphous carbon and magnetite grains grows with progen- 
itor mass. 

(iv) PISN explosions lead to the formation of grains with 
typical sizes which range between 10~^ and a few 0.1/^m; for 
all progenitors, grains radii are always smaller than 1 iim but 
the relative grain size distribution change depending on the 
thermodynamics and abundance of the relevant key species. 
As a general trend, corundum represents the smallest-size 
compound, amorphous carbon grains tend to form with the 
largest radii and silicates and magnetite grains have inter- 
mediate sizes. 

The above results have been obtained under the as- 
sumption that the heavy elements present in the ejecta at 
the onset of the explosion are uniformly mixed up to the 
outer edge of the helium core. In previous dust formation 
models (Kozasa et al. 1989; Todini & Ferrara 2001) this 
approximation was motivated by the eaxly emergence of X- 
rays and 7-rays in SN 1987A, which suggested mixing in 
the ejecta at least up to the helium core edge (Kumagai 
et al. 1988). The use of multi-dimensional hydrodynamic 
codes to model the observed light curves has made clear 
that mixing occurs on macroscopic scales through the de- 
velopment of Rayleigh- Taylor instabilities: these instabilities 
arise at the interface between elemental zones and grow non- 
linearly to produce (i) fingers of heavy elements projected 
outwards with high velocities and (ii) mixing of lighter ele- 
ments down to regions that have lower velocities (Wooden 
1997). As a result, the gas in the ejecta is mixed into regions 
which are still chemically homogeneous and which cool with 
different timescales, whereas only small clumps in the ejecta 
are microscopically mixed. Such a structure affects the pro- 
cess of dust formation, changing both the total amount of 
dust formed and the relative abundance of different solid 
compounds, as recently shown by Nozawa et al. (2004). 

A related aspect of the model which requires deeper in- 
vestigation is the assumed temperature structure of the gas 
within the ejecta. The adopted thermal evolution of the ex- 
panding ejecta is motivated by the results of numerical sim- 
ulations of the internal structure of PISN progenitors (Fryer, 
Woosley & Heger 2001; Heger & Woosley 2002) and mod- 
els the impact of the decays of radioactive elements through 
the change of the adiabatic index from 7 = 4/3 to 7 = 5/3 
assumed to occurr at ^®Co e- folding time (111.26 days after 
the explosion). A comparable thermal history is found by 
Nozawa et al. (2003) solving the radiative transfer equation 
and taking into account the energy deposition by radioactive 
elements. For the same PISN progenitor mass, however, their 
normalization is slightly higher (by a factor ~ 3) because of 
the different input stellar models (Umeda & Nomoto 2002). 

Finally, to quantify the cosmic relevance of dust forma- 
tion in the early Universe, we should restrict the analysis to 
the fraction of the newly formed dust which is able to sur- 
vive the impact of the reverse shock, following thereafter the 
fate of the surrounding metal-enriched gas. The process of 
dust sputtering dissociates the grains into their metal com- 
ponents and might have an important role in cosmic metal 
enrichment (Bianchi & Ferrara 2004). However, new theo- 
retical models seem to indicate that the post-shock temperan- 
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ture enters the regime suitable for dust condensation inside 
the oxygen layer, because of the high cooling efficiency of 
this element, but remains substantially higher in the outer 
He and H layers. If this is the case, then the reverse shock 
might lead to an increase in the total amount of dust formed 
rather than decreasing it. 



5 COSMOLOGICAL IMPLICATIONS 

Our analysis shows that if the first stars formed according 
to a top-heavy IMF and a fraction of them exploded as pair- 
creation supcrriovac, a largo amount of dust is produced in 
the early Universe. In this section, we discuss some of the 
main cosmological implications of this result, with particu- 
lar emphasis on its role in the thermodynamics of the gas 
that will be later incorporated into subsequent generations 
of objects. 

Recent numerical and semi-analytical models for the 
collapse of star- forming gas clouds in the early Universe have 
shown that because of the absence of metals and the reduced 
cooling ability of the gas, the formation of low-mass stars is 
strongly inhibited. In particular, below a critical threshold 
level of metallicity of Z^v = 10~^^^ Zq cooling and fragmen- 
tation of the gas clouds stop when the temperature reaches a 
few hundreds K (minimum temperature for II2 cooling) and 
the corresponding Jeans mass is of the order of 10^ - 10"* M© 
(Schneider et al. 2002). Gas clouds with mass comparable to 
the Jeans mass start to gravitationally collapse without fur- 
ther fragmentation, until a central protostellar core is formed 
which rapidly grows in mass through gas accretion from the 
surrounding envelope. The absence of metals and dust in 
the accretion flow and the high gas temperature favour very 
high accretion efficiencies and the resulting stars can be as 
massive as 600 Mq (Omukai & Palla 2003). 

As the gas becomes more and more enriched with heavy 
elements, the cooling rate increases because of metal (espe- 
cially C and O) line emission. More importantly, if a fraction 
of the available metals is depleted onto dust grains, dust-gas 
thermal exchages activate a new phase of cooling and frag- 
mentation which enables the formation of gas clumps with 
low-mass (Schneider et al. 2003). In particular, if the metal- 
licity of the star forming gas clouds exceeds 10~'^Zq, this 
dust-driven cooling pathway is irrelevant because cooling via 
metal-line emission by itself is able to fragment the gas down 
to characteristic Jeans masses in the range 10"^ - IMq. At 
the same time, if the metallicity of the gas is below W~^Zq, 
even if all the available metals are assumed to be depleted 
onto dust grains, the resulting cooling efficiency is too low to 
activate fragmentation and the resulting Jeans masses arc as 
large as in the metal-free case (10=^- 10'' Mq). Thus, we can 
conclude that the presence of dust is crucial to determine 
the final mass of stars forming out of gas clouds with metal- 
licitics in the critical range Za — 10"^^^^©. As a reference 
value, for a metallicity oi Z = 10^'''^ Zq, if 20% of the met- 
als are depleted onto dust grains (/dep = 0.2) the resulting 
mass of the collapsing clouds is reduced to 10~^ — 10~^Mo 
and can lead to low-mass stars. This is fully consistent with 
the predicted dust mass formed from PISN progenitors (see 
Fig. 6). 

Also, we can not neglect that our models predict the 
formation of a significant amount of CO molecules, which 



can contribute to gas cooling. This aspect needs to be in- 
vestigated further though we expect that the presence of 
CO molecules will be complementary to that of dust in the 
critical range of metallicities. 

Therefore, in the emerging picture of galaxy evolution, 
the first episodes of star formation were characterized by 
very massive stars forming according to a top-heavy IMF. 
It is only when the metals and dust ejected by the first 
pair-creation supernovae are able to pollute a substantial 
fraction of the IGM, that an overall transition to a normal 
IMF forming stars with masses comparable to those that we 
presently observe in the nearby Universe occurs. The epoch 
of this transition depends crucially on the process of metal 
enrichment through the so-called chemical feedback: within 
metal-enriched regions, if the metallicity lies within the crit- 
ical range, the presence of dust becomes critical and can no 
longer be neglected in the thermodynamics of the star form- 
ing gas. It is very likely that the transition will not occur at 
a single redshift because of the highly inhomogeneous pro- 
cess of metal enrichment (Scannapieco, Schneider & Ferrara 
2003) and that there will be epochs when the two modes 
of star formation will be coeval in different regions of the 
Universe. 

This in turns might bo very important for the rcion- 
ization history of the IGM. Indeed, very massive metal-free 
stars are powerful sources of ionizing photons and an early 
epoch of star formation with a top-heavy IMF can easily 
match the required high optical depth to electron scattering 
measured by WMAP (Kogut et al. 2003; Spergel et al. 2003; 
for a general discussion see Ciaxdi, Ferraxa & White 2003). 
However, an early epoch of reionization is difficult to recon- 
cile with the observed Gunn-Peterson effect in the spectra 
of z > 6 quasars (Fan et al. 2002) which implies a mass 
averaged neutral fraction of ~ 1%. These two observational 
constraints seem to indicate that the reionization history of 
the Universe might have been more complex than previously 
thought, with probably two distinct epochs of reionization 
separated by a relatively small redshift interval in which the 
Universe might have recombined again. This intermediate 
epoch of recombination could be the result of a decrease 
in the ionizing power of luminous sources, probably due to 
their metal-enrichment and consequent IMF transition. 

An early epoch of reionization, as required by WMAP 
data, poses another critical issue: after reionization, the tem- 
perature of the IGM starts to decline as a consequence of cos- 
mic expansion. By the time it reaches the observable range 
of redshifts z < 6 the temperature of the IGM might be 
too low to match the observed thermal history (Schayc et 
al. 2000). The presence of a significant amount of dust syn- 
thetized by the first very massive supernovae might be ex- 
tremely important to raise the IGM temperature through 
dust photoelectric heating. To estimate the amount of dust 
in the IGM required to match the observed thermal history, 
it is necessary to make specific assumptions about the UV 
background radiation, the reionization history of the IGM 
and the grain composition and size distribution (Inoue & 
Kamaya 2003). Furthermore, the properties of dust grains 
in the IGM may differ from those directly predicted by dust 
formation models as a consequence of specific selection rules 
in the transfer of grains from the host galaxies to the IGM 
(see Bianchi & Ferrara 2004, who find that only grains larger 
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Figure 8. The number of PISN required to pollute a host proto- 
galaxy witli a dust-to-gas ratio equal to 5 % of the galactic value 
is shown as a function redshift (solid line). In the same panel, we 
show the total mass of the host halo which correspond to a virial 
temperature of 1000 K (dashed line). 



than ~ 0.06/im arc preferentially ejected in the IGM). These 
complications are beyond the scope of the present analysis. 

Finally, it is well known that the presence of dust at 
high redshifts offers an alternative formation channel for 
molecular hydrogen, the dominant coolant in the early Uni- 
verse, which, in the absence of dust, can form only from 
the gas phase. In small protogalaxies, the H2 formation rate 
on grain surface becomes dominant with respect to the for- 
mation rate from the gas phase, when the dust-to-gas ratio 
exceeds roughly 5% of the galactic value (Todini & Ferrara 
2001). From this order of magnitude estimate, we can derive 
the number of pair-creation supernovae required to enrich a 
small protogalaxy to this level. The results are shown in 
Fig. 8, where we have assumed the host galaxy to corre- 
spond to a halo with virial temperature of 1000 K and we 
have considered the dust mass synthetized by a 195 Mq 
progenitor (see Fig. 4) as representative of dust formation 
efficiencies in PISN. This simple estimate shows that in the 
redshift range relevant to these small protogalactic systems 
{z > 10 — 15), one PISN is required to pollute the gas with 
enough dust that H2 formation on grain surface starts to 
be important. This, in turns, might have very important 
consequences for the star formation activity at high redshift 
(Hirashita & Ferrara 2002). 
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